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Abstract. We present a detailed optical spectroscopic and B, V , I, Ha photometric study of the metal-deficient 
cometary blue compact dwarf (BCD) galaxy SBS 1415+437. We derive an oxygen abundance 12 + log(0/H) = 
7.61±0.01 and 7.62±0.03 (Z = Z Q /20)* in the two brightest H II regions, among the lowest in BCDs. The helium 
mass fractions in these regions are Y — 0.246±0.003 and 0.243±0.010. Four techniques based on the equivalent 
widths of the hydrogen emission and absorption lines, the spectral energy distribution and the colours of the galaxy 
are used to put constraints on the age of the stellar population in the low-surface-brightness (LSB) component of 
the galaxy, assuming two limiting cases of star formation (SF), the case of an instantaneous burst and that of a 
continuous SF with a constant or a variable star formation rate (SFR). The spectroscopic and photometric data 
for different regions of the LSB component are well reproduced by a young stellar population with an age t < 250 
Myr, assuming a small extinction in the range Av — - 0.6 mag. Assuming no extinction, we find that the upper 
limit for the mass of the old stellar population, formed between 2.5 Gyr and 10 Gyr, is not greater than ~ (1/20 
- 1) of that of the stellar population formed during the last ~ 250 Myr. Depending on the region considered, this 
also implies that the SFR in the most recent SF period must be 20 to 1000 times greater than the SFR at ages 
i; 2.5 Gyr. We compare the photometric and spectroscopic properties of SBS 1415+437 with those of a sample 
of 26 low-metallicity dwarf irregular and BCD galaxies. We show that there is a clear trend for the stellar LSB 
component of lower-metallicity galaxies to be bluer. This trend cannot be explained only by metallicity effects. 
There must be also a change in the age of the stellar populations. The most metal-deficient galaxies have also 
smaller luminosity-weighted ages. 
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1. Introduction 

Since its discovery as a metal-deficient blue compact dwarf 
(BCD) galaxy (Thuan, Izotov & Lipovetskv HM5|) . SBS 
1415+437 (= CG 389) has been considered as a probable 
nearby young dwarf galaxy. Situated at a distance D = 
11 A Mpc it was classified by Thuan, Izotov & Foltz l|1999[) 
as a cometary BCD with a very bright supergiant H II 
region at the SW tip of the galaxy. 

From 4m Kitt Peak National Observatory (KPNO) 
telescope spectra, Thuan et al. I|1995fl first derived an 
oxygen abundance of 12 + log(0/H) = 7.51+0.01 in 
SBS 1415+437 placing the galaxy among the most metal- 
deficient BCDs known. Later, Izotov & Thuan (1998 
1999) derived from the same spectrum 12 + log(0/H) = 

Send offprint requests to: N.G. Guseva, guseva@mao.kiev.ua 
* 12+log(O/H) = 8.92 (Anders & Grevesse 115551 . 



7.59+0.01 using five-level atom models for abundance de- 
termination instead of the three- level atom model used by 
Thuan et al. iflMSjl . Thuan et al. l(T^fjjl . using Multiple 
Mirror Telescope (MMT) and Hubble Space Telescope 
(HST) FOS observations, derived 12 + log(0/H) = 
7.60+0.01 and 7.54+0.14, respectively. The high bright- 
ness and low metallicity of the H n region in SBS 
1415+437 make this galaxy one of the best objects for 
helium abundance determination. Izotov & Thuan (1998) 
and Thuan et al. 1)1999(1 derived respectively a helium mass 
fraction Y = 0.244+0.002 and 0.246+0.004 for it, close to 
the primordial helium mass fraction of Y p = 0.245+0.002 
by Izotov et al. fllMty . 

Thuan et al. ( 1993} have discussed the evolutionary 
status of SBS 1415+437, using ground-based MMT spec- 
troscopic and iTST/WFPC2 photometric data. Based on 
the (V — I) vs I colour-magnitude diagrams (CMD) and 
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spectral energy distributions (SED) in the optical range, 
they concluded that SBS 1415+437 is a truly young galaxy 
that did not start to form stars until ~ 100 Myr ago. 
However, the V and I images used by Thuan et al. (1999) 
were not deep enough for the detection of old red giant 
branch (RGB) stars in the CMD. Furthermore, they con- 
sidered an instantaneous burst model which gives only a 
lower limit to the age of the stellar population in SBS 
1415+437. 

In this paper we combine new spectroscopic and photo- 
metric data with previous observations to derive elemental 
abundances and to better constrain the age of the stellar 
population in SBS 1415+437. For the latter task we use 
four different techniques of age determination and con- 
sider different star formation (SF) histories. The paper 
is organized as follows. In Sect. [21 we describe the obser- 
vations and data reduction. The photometric properties 
of SBS 1415+437 are described in Sect. El We derive in 
Sect. 01 the chemical abundances in the two brightest H II 
regions. In Sect. we discuss the properties of the stel- 
lar populations in SBS 1415+437 and compare them with 
those in other low-met allicity dwarf galaxies. Finally, Sect. 
El summarises the main conclusions of this study. 

2. Observations and data reduction 

2.1. Photometry 

Narrow-band images of SBS 1415+437 in the Ha line at 
A6563A through a passband with a full width at half max- 
imum (FWHM) of 74A, and in the adjacent continuum at 
A6477A through a passband with FWHM = 72A were 
obtained with the Kitt Peak 1 2.1m telescope on April 
22, 1999 during a photometric night. The telescope was 
equipped with a Tektronix 1024 x 1024 CCD detector 
operating at a gain of 3e~ ADU -1 , giving an instrumen- 
tal scale of 0'.'305 pixel -1 and field of view of 5'. The total 
exposures of 50 min in the Ha line and 40 min in the adja- 
cent continuum bluewards of Ha were split up into 5 and 4 
subexposures, slightly offset with respect to each other for 
removal of cosmic particle hits and bad pixels. The point 
spread function has a FWHM of 2'.'2. Bias and flat-field 
frames were obtained during the same night. The standard 
stars Feige 34 and HZ 44 were observed in both filters dur- 
ing the same night at several airmasses for absolute flux 
calibration. 

Another broad-band B image (15 min) of SBS 
1415+437 was obtained on March 9, 1997 under pho- 
tometric conditions, with the CAFOS focal reducer at- 
tached to the 2.2m telescope of the German-Spanish 
Astronomical Center, Calar Alto 2 , Spain. CAFOS was 

1 Kitt Peak National Observatory (KPNO) is operated by 
the Association of Universities for Research in Astronomy 
(AURA), Inc., under cooperative agreement with the National 
Science Foundation (NSF). 

2 German-Spanish Astronomical Center, Calar Alto, oper- 
ated by the Max-Planck-Institute for Astronomy, Heidelberg, 
jointly with the Spanish National Commission for Astronomy. 
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Fig. 1. HST V image of SBS 1415+437. The long-slit 
positions during the two observations with the MMT and 
4m Mayall telescope are labeled "1" and "2" respectively. 
North is up and east is to the left. Regions el - e3 with 
hydrogen emission lines in the spectra and regions al - a4 
with hydrogen absorption lines are labeled. 

equipped with a SITe 2048 x 2048 CCD operating at a 
gain of 2.3 e~ ADU^ 1 , with a read-out noise of < 3 counts 
(rms). With a focal ratio of f/4.4, the instrumental scale 
was 0'.'53 pixel^ 1 and field of view ~ 15'. The seeing during 
the observations was 3'.'1 (FWHM). Standard stars from 
Christian et al. (1985) were observed at several airmasses 
for calibration. 

Standard reduction steps, including bias subtraction, 
flat-field correction, removal of cosmic particle hits and 
absolute flux calibration were carried out using IRAF 3 
and MIDAS 4 . 

The ground-based photometric data were supple- 
mented by the HST/WFPC2 V (F569W) and I (F791W) 
images described by Thuan et al. ([1999 ) . 

2.2. Spectroscopy 

Spectroscopic observations were carried out on June 18, 
1999, at the Kitt Peak 4m Mayall telescope with the 
Ritchey- Chretien spectrograph and a T2KB 2048 x 2048 

3 IRAF is the Image Reduction and Analysis Facility dis- 
tributed by the National Optical Astronomy Observatory, 
which is operated by the AURA under cooperative agreement 
with the NSF. 

4 Munich Image Data Analysis System, provided by the 
European Southern Observatory (ESO). 
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Fig. 2. Continuum-subtracted Ha image of SBS 
1415+437 with superposed Ha continuum isophotes. 
Regions with nebular emission are labeled el, e2 and e4. 
The faint region e3 with H/3 and Ha emission lines in its 
spectrum is not seen in the Ha image. 



CCD detector. The 2" x 300" slit was centered on the 
brightest H n region el (slit 2 in Fig.^) with position angle 
P.A. = 48° so as to include the second brightest H n region 
e2 to the SW of region el. We used the KPC-10A grating 
in first order and a GG 375 order separation filter. The 
spatial scale along the slit was 0'/69 pixel -1 and the spec- 
tral resolution ~7 A (FWHM). The spectra were obtained 
at an airmass 1.27. The total exposure time of 60 minutes 
was broken up into 3 subexposures. No correction for at- 
mospheric refraction was made because of the small air- 
mass during the observations. Two Kitt Peak spectropho- 
tomctric standard stars were observed for flux calibration. 
For wavelength calibration, He-Ne-Ar comparison spectra 
were obtained after each exposure. 

The data reduction was performed with the IRAF 
software package. This includes bias-subtraction, flat- 
field correction, cosmic-ray removal, wavelength calibra- 
tion, night sky background subtraction, correction for at- 
mospheric extinction and absolute flux calibration of the 
two-dimensional spectrum. 

For abundance determination, one-dimensional spec- 
tra of regions el and e2 were extracted within apertures 
of 2" x 4'.'6 and 2" x 4'.'0, respectively. In addition, we 
extracted spectra of the low-surface-brightness (LSB) re- 
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Fig. 3. (a) Surface brightness profiles (SBPs) of SBS 
1415+437 in V and I (filled and open circles, respectively), 
derived from H5T/WFPC2 data. A linear fit to the V 
SBP in the radius interval 4"< R* <T3" is shown by the 
upper solid-grey line. Both SBPs reveal a significant slope 
change for R* ^13". Filled squares show the ground-based 
B SBP, shifted vertically by 0.5 mag. In addition to the 
exponential intensity slope at radii 4"< R* ^13", this 
profile reveals for R* >16" an outer exponential regime 
with an a ~0.27 kpc (lower solid-grey line). The residu- 
als after subtraction of the outer exponential distribution 
from the B SBP are shown by crosses, (b) (V — I) colour 
profile of SBS 1415+437, computed from subtraction of 
the SBPs in (a). A strong colour gradient j + = 2.5 mag 
kpc -1 is derived for radii smaller than the effective radius 
r e fl (inner fit). The (V — I) colour of 0.4 - 0.5 mag for 
larger radii is nearly constant with a 7+ <0.2 mag kpc" 1 . 
Comparison of the B and V SBPs implies a B — V color 
-0.2 mag in the radius range 13" < R* < 16". 



gions a3 and a4 with strong hydrogen Balmer absorption 
lines, and of region e3 with Ha and H/3 in emission (Fig. 

ID 

We also used the two-dimensional MMT spectrum ob- 
tained by Thuan et al. I|1999f) with the slit oriented at 
P.A. = 22° (slit 1 in Fig. [I}. We extracted one-dimensional 
spectra of the LSB regions al and al within apertures of 
T/5 x 3"4 and l'/5 x 13"2, respectively (Fig. UJ). These 
spectra show strong hydrogen Balmer absorption lines. 
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Fig. 4. The KPNO 4m telescope spectra of the brightest H II regions el and e2 with the identified emission lines. The 
lower spectra in (a) and (b) are the observed spectra downscaled by factors of 50 and 30, respectively. 



The selected LSB regions are listed in Tables and 0] 
with their positions and spatial extents. Origins are set 
at the center of the brightest region el (Fig. [2J . The 
spectra of the LSB regions are used to study the stellar 
populations and to constrain the age of the oldest stars 
which contribute to the light of these regions. 

3. Photometric properties 

It is seen from the continuum-subtracted Ha images (Fig. 
2) that star-forming activity in SBS 1415+437 is primarily 
occurring in regions el and e2, with some additional faint 
Ha emission present in region e4. However, the available 
narrow-band data are not deep enough for tracing faint 



Ha emission in other regions of the LSB component, such 
as in regions al, a3 and e3 where Ha and H/? have been 
detected spectroscopically (Table 0}. The latter fact sug- 
gests that some low-level SF may be present at various 
locations within the LSB component. 

The photometric properties of the stellar LSB compo- 
nent of SBS 1415+437 were first investigated by Thuan 
et al. (1999) using HST/WFPC2 V and I images. These 
authors fit the surface brightness profiles (SBPs) of SBS 
1415+437 with an exponential distribution in the radius 
range 4"< R* £13" with a scale length a = 5'.'4 (w0.3 
kpc). However, their study was limited by the fact that 
the HST images do not include the outermost NE part of 
the LSB component. Furthermore, the SBPs in Thuan et 
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al. (1999) reach only a surface brightness level /i ~ 24.5 
mag arcsec -2 . It is known, however, that the star- forming 
component may contribute to the optical BCD emission 
to fainter surface brightness levels (see e.g. Papaderos et 
al. 2002 and references therein). A comparatively young 
stellar population has been observed in the inner part of 
the LSB component of other cometary BCDs, several hun- 
dred pc away from the brightest H n region (Noeske et al. 
2000; Guseva et al. 2001; Fricke et al. 2001). Therefore 
deeper images are needed to study the outer parts of the 
LSB component in SBS 1415+437. 

Note that the SBPs by Thuan et al. (1999) show in 
the outermost part (13" < R* < 16", or 24 < \x v < 24.5) 
a steeper exponential intensity decrease than the one ob- 
served at intermediate intensity levels. This slope change, 
not discussed in Thuan et al. (1999), is found indepen- 
dently by us in the H5T/WFPC2 SBPs derived with the 
method iv of Papaderos et al. (2002) and ellipse fitting 
to the visible part of the LSB component. While the lat- 
ter method extends surface photometry out to larger radii 
(R* ~ 22"), it is subject to large uncertainties because 
the NE part of the LSB component with [i ^24.5 V mag 
arcsec" 2 (R* > 16") lies outside the JSST/WFPC2 field ° 
of view. Additionally, the outermost LSB isophotes show ^ 
considerable deviations from ellipticity. 

To study the surface brightness distribution at large 
radii we use the ground-based Calar Alto B image. Despite \ 
the poor spatial resolution this image allows us to study £ 
the entire LSB component out to its Holmberg radius. The 
change in the exponential slope for 13" Ss R* Ss 16" is con- 
firmed from the ground-based B SBP. At large radii, how- o 
ever, this SBP reveals a flatter, outer exponential regime 
with a scale length a fairly comparable to that previously 
obtained at intermediate intensity levels from HST data 
(upper thick-grey line in Fig. 3a). From fitting an expo- 
nential model to the B SBP for R* >16" (lower thick- 
grey line in Fig. 3a) we obtain a central surface brightness 
A*e,o = 21.37 B mag arcsec -2 and a scale length a = 0.27 
kpc. 

Note, however, that the inner exponential profile stud- 
ied by Thuan et al. (1999) is J>0.3 mag brighter than the 
outer one, which suggests that more than 1/4 of the emis- 
sion associated with this profile originates from the part 
of the LSB component between regions el and al. 

The present data provide no compelling evidence for a 
large age difference between the stellar population which 
dominates within the inner exponential regime discussed 
in Thuan et al. (1999) and that responsible for the outer- 
most LSB emission (i.e. for R* > 16"). The {V - 1) profile 



reveals a strong colour gradient (7 + =2.5 mag kpc , inner 
solid-grey line in Fig. Et>) within the inner 5", or roughly 
the V band effective radius (r off = 5'.'6) of SBS 1415+437. 
At larger radii, however, linear fits to the (V — /) pro- 
file yield, depending on whether they are error-weighted 
or not, a gradient not exceeding 0.1 and 0.2 mag kpc -1 , 
respectively. Comparable values are also found from sub- 
traction of the exponential fits in V and I (Fig. UJi) in the 



radius range r e g < R* < 13" (Fig. |2t>, solid-grey line at 
intermediate radii). 

4. Chemical abundances 

In this section we derive the elemental abundances of re- 
gions el and e2 using the Kitt Peak 4m telescope observa- 
tions. Their spectra with strong emission lines are shown 
in Fig. E| 

The observed (F(X)) and extinction-corrected (/(A)) 
emission line fluxes relative to the H/3 emission line fluxes, 
their equivalent widths EW, the extinction coefficients 
C(H/3), the observed fluxes of the H/3 emission line, and 
the equivalent widths of the hydrogen absorption lines for 
regions el and e2 are shown in Tabled Despite the dif- 
ferences in aperture (2" x 4" 6 for the Kitt Peak 4m data, 
1"5 x 5" and 1"5 x 0'.'6 for the MMT data from Thuan 
et al. 1999), the relative fluxes of the emission lines for 
region el are in agreement within the errors with those 
derived by Thuan et al. I|1999(l . 
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Fig. 5. Comparison of the elemental abundance ratios, ob- 
tained for the brightest H n regions el (large squares) and 
e2 (small squares) with data for other BCDs (open circles). 

The physical conditions and heavy element abun- 
dances in regions el and e2 were derived following Izotov 

mm 



et al. l(raMllT997a|l and Thuan et al. (|T995l) . The electron 
temperatures T e (0 in), T (S in), T e (0 il) for the high-, 
intermediate- and low-ionization regions respectively, the 
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Table 1. Observed (F(X)) and extinction-corrected (I(X)) fluxes and equivalent widths (EW) of emission lines in the 
H II regions el and e2. 



Ao(A) Ion 


region el 




region e2 




F(A)/F(H/3) J(A)/J(H/J) 


EW(A) 


F{\)/F{R(3) I{\)/I(H0) 


EW(A) 



3727 [O n] 







.968 


±0 


.004 


0.949 


±0.004 


88 


.2 


±0 


.3 


1. 


.099 


±0. 


.016 


1.085 


±0.016 


79. 


.2 


±0. 


.8 


3750 H12 







.024 


±0 


.001 


0.060 


±0.004 


2 


.2 


±0 


.1 


0. 


.007 


±0. 


.004 


0.031 


±0.021 


0. 


.5 


±0. 


.3 


3771 Hll 







.029 


±0 


.001 


0.064 


±0.003 


2 


.7 


±0 


.1 


0. 


.030 


±0. 


.004 


0.053 


±0.009 


2 


.3 


±0. 


.3 


3798 H10 







.038 


±0 


.001 


0.073 


±0.003 


3 


.5 


±0 


.1 


0. 


.045 


±0 


.005 


0.067 


±0.009 


3. 


.4 


±0. 


.4 


3835 H9 







.057 


±0 


.001 


0.091 


±0.003 


5 


.3 


±0 


.1 


0. 


.055 


±0 


.005 


0.078 


±0.008 


4. 


.2 


±0. 


.4 


3868 [Ne in] 







.262 


±0 


.002 


0.257 


±0.002 


25 


.2 


±0 


.2 


0. 


.241 


±0 


.007 


0.238 


±0.007 


18. 


.4 


±0. 


.5 


3889 H8 + He i 





.166 


±0 


.002 


0.196 


±0.002 


16 


.4 


±0 


.2 


0. 


.151 


±0. 


.006 


0.172 


±0.008 


11 


.6 


±0. 


.4 


3968 [Ne in] 


+ H7 





.227 


±0 


.002 


0.255 


±0.003 


23 


.0 


±0 


.2 


0. 


.217 


±0. 


.007 


0.236 


±0.008 


17 


.2 


±0. 


.5 


4026 He I 







.013 


±0 


.001 


0.013 


±0.001 


1 


.3 


±0 


.1 






















4068 [S n] 







.008 


±0 


.001 


0.008 


±0.001 





.8 


±0 


.1 






















4101 H5 







.245 


±0 


.002 


0.269 


±0.002 


27 


.6 


±0 


.2 


0. 


.236 


±0. 


.007 


0.253 


±0.008 


19. 


.9 


±0. 


.5 


4340 H 7 







.452 


±0 


.002 


0.469 


±0.003 


56 


.9 


±0 


.3 


0. 


.446 


±0. 


.009 


0.458 


±0.010 


42 


.9 


±0. 


.7 


4363 [O in] 







.085 


±0 


.001 


0.084 


±0.001 


10 


.8 


±0 


.2 


0. 


.071 


±0. 


.005 


0.070 


±0.005 


6. 


.9 


±0. 


.5 


4388 He I 







.004 


±0 


.001 


0.004 


±0.001 





.5 


±0 


.2 






















4471 He i 







.034 


±0 


.001 


0.033 


±0.001 


4 


.5 


±0 


.2 


0. 


.034 


±0. 


.004 


0.034 


±0.004 


3. 


.6 


±0. 


.4 


4658 [Fe in] 







.006 


±0 


.001 


0.006 


±0.001 


1 


.0 


±0 


.1 






















4686 He n 







.023 


±0 


.001 


0.023 


±0.001 


3 


.5 


±0 


.2 






















4713 [Ar iv] 


+ He I 





.011 


±0 


.001 


0.011 


±0.001 


1 


.7 


±0 


.2 






















4740 [Ar iv] 







.005 


±0 


.001 


0.005 


±0.001 





.8 


±0 


.2 






















4861 R/3 




1 


.000 


±0 


.004 


1.000 


±0.004 


166 


.2 


±0 


.5 


1. 


.000 


±0. 


.014 


1.000 


±0.014 


133 


.8 


±1. 


.3 


4922 He I 







.008 


±0 


.001 


0.008 


±0.001 


1. 


.4 


±0 


.2 


0. 


.012 


±0. 


.004 


0.012 


±0.004 


1. 


.7 


±0. 


.5 


4959 [O in] 




1 


.203 


±0 


.005 


1.179 


±0.005 


214. 


.2 


±0 


.6 


1. 


.082 


±0. 


.015 


1.068 


±0.015 


156 


.0 


±1. 


.5 


5007 [O in] 




3 


.613 


±0 


.012 


3.542 


±0.012 


663 


.0 


±0 


.9 


3. 


.210 


±0. 


.036 


3.169 


±0.036 


473. 


.6 


±2. 


.4 


5199 [N i] 







.005 


±0 


.001 


0.005 


±0.001 


1 


.0 


±0 


.2 






















5271 [Fe in] 







.003 


±0 


.001 


0.003 


±0.001 





.7 


±0 


.2 






















5518 [CI in] 







.003 


±0 


.001 


0.003 


±0.001 





.8 


±0 


.2 






















5538 [CI in] 







.002 


±0 


.001 


0.002 


±0.001 





.4 


±0 


.2 






















5876 He I 







.100 


±0 


.001 


0.098 


±0.001 


28 


.8 


±0 


.3 


0. 


.107 


±0 


.004 


0.106 


±0.004 


26 


.2 


±1. 


.0 


6300 [O i] 







.022 


±0 


.001 


0.022 


±0.001 


8 


.0 


±0 


.3 


0. 


.032 


±0. 


.003 


0.032 


±0.003 


9. 


.8 


±0. 


.9 


6312 [S in] 







.015 


±0 


.001 


0.015 


±0.001 


5 


.3 


±0 


.2 


0. 


.012 


±0. 


.002 


0.011 


±0.002 


3 


.5 


±0. 


.7 


6364 [O i] 







.008 


±0 


.001 


0.008 


±0.001 


2 


.7 


±0 


.2 


0. 


Oil 


±0. 


.002 


0.011 


±0.002 


3 


.4 


±0. 


.7 


6563 Ha 




2 


.714 


±0 


.009 


2.670 


±0.010 


997. 


.9 


±1 


.7 


2. 


.810 


±0 


.031 


2.780 


±0.034 


871. 


.7 


±4. 


.5 


6584 [N n] 







.033 


±0 


.001 


0.032 


±0.001 


11 


.1 


±0 


.3 


0. 


.042 


±0. 


.003 


0.042 


±0.003 


13. 


.2 


±1. 


.0 


6678 He I 







.029 


±0 


.001 


0.028 


±0.001 


11 


.1 


±0 


.3 


0. 


.029 


±0 


.003 


0.028 


±0.003 


9 


.2 


±0. 


.8 


6717 [S n] 







.090 


±0 


.001 


0.088 


±0.001 


34 


.1 


±0 


.4 


0. 


.107 


±0. 


.004 


0.105 


±0.004 


34 


.3 


±1. 


.2 


6731 [S ii] 







.066 


±0 


.001 


0.065 


±0.001 


25 


.2 


±0 


.4 


0. 


.080 


±0. 


.004 


0.079 


±0.004 


25 


.9 


±1. 


.1 


7065 He I 







.027 


±0 


.001 


0.026 


±0.001 


11 


.8 


±0 


.3 


0. 


.027 


±0. 


.002 


0.027 


±0.002 


8. 


.4 


±0. 


.9 


7136 [Ar m] 







.043 


±0 


.001 


0.042 


±0.001 


20 


.6 


±0 


.4 


0. 


.048 


±0 


.003 


0.048 


±0.003 


21. 


.0 


±1. 


.3 


C(H/3) dex 












0.000±0.004 


















0.000±0.014 










F(H/3) a 












4.54±0.01 


















0.71±0.01 










£W(abs) A 












3.4±0.1 


















1.8±0.3 











Mn units 10 14 erg s 'cm 2 . 



electron number densities N e (S Ii), ionization correction 
factors (ICF), and ionic and total heavy element abun- 
dances are shown in Tabled for both regions. The oxygen 
abundance 12 + log(0/H) = 7.61 ± 0.01 (Z Q /20) and 
heavy element abundance ratios for region el are in good 
agreement with those derived by Thuan et al. I|1999|l ■ The 
oxygen abundance 12 + log(0/H) = 7.62 ± 0.03 and heavy 
element abundance ratios in region e2 are consistent with 
those for region el within the errors. 

In Figure |3 we compare the heavy element abundance 
ratios in the two brightest regions of SBS 1415±437 with 
data for a sample of low-metallicity BCDs. The Ne/O, 



S/O, Ar/O and [O/Fe] abundance ratios for the compar- 
ison sample are taken from Izotov & Thuan (1999), while 
the Cl/O abundance ratios are collected from Izotov & 
Thuan (1998) and Izotov et al. (1997a). The heavy element 
abundance ratios for regions el (large squares) and e2 
(small squares) are in good agreement with those for other 
BCDs. Note, that the Cl/O ratio does not show any sig- 
nificant increase with increasing oxygen abundance. This 
conclusion is strengthened by the observations of Esteban 
et al. (1998, 1999a, 1999b) who derived log(Cl/0) in the 
range from -3.28 to -3.47 for high-metallicity H II regions 
in Orion, M17 and M8 with 12+log(0/H) = 8.60, 8.50 
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Table 2. Element abundances in regions el and e2. 



Value 



region el 



region el 



4000 4500 
Wavelength (A) 



Fig. 6. The blue part of the MMT spectrum of region al 
with labeled emission and absorption lines. 



and 8.60, respectively. For comparison, log(Cl/0) = -3.37 
is derived for the Sun (Anders & GrevesseQ989). 

The high brightness of regions el and e2 allows for 
a reliable determination of the 4 He abundance. Nine He 
I emission lines are detected in the spectrum of region 
el (Table P). Two of them, He I A3889 and A4713, are 
blended with other lines. Six He I lines are detected in 
region e2. The five brightest He I A3889, A4471, A5876, 
A6678, A7065 emission lines are used to correct their fluxes 
for collisional and fluorescent enhancement. This is done 
by minimizing the deviations of the corrected He I line 
flux ratios from the recombination ratios, through vary- 
ing the electron number density in the He + zone and the 
optical depth in the He I A3889 emission line. The flux of 
this line was preliminarily corrected for the contribution of 
the H I A3889 emission line, according to prescriptions of 
Izotov et al. (1994, 1997a). Helium abundances He + /H + , 
derived from the corrected He I A4471, A5876, A6678 line 
fluxes and their weighted mean are listed in Table [21 The 
abundance He +2 /H + is added to He + /H + for region el, as 
He II A4686 is present in its spectrum. Note the lower He 
abundance derived from the He I A4471 flux which is most 
likely due to significant underlying stellar He I A4471 ab- 
sorption. The effect of underlying absorption for the other 
He I emission lines used in the He abundance determina- 
tion is much smaller, as they have much larger equivalent 
widths compared to the He I A4471 emission line (Table 
HJ). The mean 4 He mass fractions Y = 0.246±0.003 and 
0.243±0.010 in regions el and e2 (Table 0) are consistent 
with the values derived for SBS 1415+437 by Izotov & 
Thuan fiMty and Thuan et al. <fT59"9"j) . They are also con- 
sistent with the primordial 4 He mass fraction Y p — 0.244 
± 0.002, derived by extrapolating the Y vs O/H linear 
regression to O/H = (Izotov & Thuan IT^5|l . or to Y p = 
0.245 ± 0.002 derived from spectroscopic observations of 
the two most metal-deficient BCDs known, I Zw 18 and 
SBS 0335-052 (Izotov et a,l. lTM?jl . 



T c (0 III)(K) 

T c (0 II) (K) 

T e (S m)(K) 
iV c (S n)(cm- 3 ) 

O+/H+(xl0 5 ) 
O +2 /H+(xl0 5 ) 
O +3 /H+(xl0 6 ) 
O/H(xl0 5 ) 
12 + log(0/H) 

N+/H+(xl0 7 ) 

ICF(N) a 

log(N/0) 

Ne +2 /H+(xl0 5 ) 

ICF(Ne) a 

log(Ne/0) 

S+/H+(xl0 7 ) 
S +2 /H+(xl0 7 ) 
ICF(S) a 
log(S/0) 

Cl ++ /H+(xl0 8 ) 

ICF(Cl) a 

log(Cl/0) 

Ar+2/H+(xl0 7 ) 
Ar+ 3 /H+(xl0 7 ) 
ICF(Ar) a 
log(Ar/0) 

Fe ++ /H+(xl0 7 ) 
ICF(Fe) a 
log(Fe/0) 
[0/Fe] b 

He+/H+(A4471) 

He+/H+(A5876) 

He+/H+(A6678) 

He+/H+ 

(weighted mean) 

He +2 /H+(A4686) 

He/H 

Y 



16490+140 
14430+110 
15380+110 
60+30 

0.930+0.020 
3.071+0.063 
0.994+0.057 
4.100+0.067 
7.61+0.01 

2.564+0.056 
4.41 

-1.560+0.017 

0.480+0.011 
1.34 

-0.806+0.014 

1.618+0.025 
6.980+0.350 
1.50 

-1.503+0.016 

1.017+0.223 
1.57 

-3.409+0.076 

1.449+0.032 
0.608+0.119 
1.05 

-2.280+0.026 

1.320+0.200 
5.51 

-1.751+0.016 
0.331+0.029 

0.071+0.002 
0.081+0.001 
0.082+0.002 

0.080+0.001 
0.002+0.000 
0.082+0.001 
0.246+0.003 



15910+540 
14170+460 
14900+450 
90+80 

1.127+0.101 
2.997+0.256 

4.123+0.275 
7.62+0.03 

3.460+0.311 
3.66 

-1.513+0.069 

0.491+0.045 
1.38 

-0.786+0.057 

2.028+0.118 
5.942+1.335 
1.39 

-1.571+0.068 



1.743+0.130 
1.64 

-2.158+0.038 



0.070+0.009 
0.083+0.004 
0.081+0.008 

0.080+0.003 

0.080+0.003 
0.243+0.010 



a ICF is the ionization correction factor. 
b [0/Fe] = log (O/Fe) - log (O/Fe) . 



5. Age of the stellar population in the LSB 
regions 

We consider next the spectroscopic and photometric prop- 
erties of the LSB regions labeled al, a2 (slit position 1 in 
Fig. flj and a3, e3, a4 (slit position 2), to constrain the 
age of the stellar populations contributing to the light in 
those regions. Ha and H/3 emission lines are present in re- 
gions al, a3 and e3 while H7 and H5 absorption lines are 
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detected in the spectra of all regions except for region e3. 
This allows us to derive the age of the stellar population 
using four methods, based on: (1) the time evolution of 
equivalent widths (EW) of hydrogen emission lines, (2) 
the time evolution of EWs of hydrogen absorption lines, 
(3) the comparison of the observed and theoretical spec- 
tral energy distributions, and (4) the broad-band colours. 
The requirement of consistency of the ages determined 
from the reddening-insensitive methods 1 and 2 and from 
the reddening-sensitive methods 3 and 4 allows to simul- 
taneously derive the extinction coefficient and constrain 
the SF history (Guseva et al. I2UUT1 !3UU5al I2003b|> . 

We measured the fluxes and equivalent widths of the 
Ha and H/3 emission lines and the H7 and US absorption 
lines in the spectra of the LSB regions, and list them in 
Tables H and H Because the H/3 emission line is narrower 
than the absorption line in these regions and does not fill 
the absorption component, its flux was measured using 
the continuum level at the bottom of the absorption line. 
This level was chosen by visually interpolating from the 
absorption line wings to the center of the line. 

The extinction coefficient C(H/3) in those regions is 
derived from the Ha/H/3 flux ratio. We adopt the theo- 
retical recombination Ha/H/3 flux ratio of 2.8, which is 
typical for hot low-metallicity H 11 regions. No correction 
for the absorption line equivalent widths has been made. 
The extinction coefficients C(Hf3) are shown in Table 

Hydrogen absorption lines are seen in the spectra of 
all regions labeled in Fig. except for the brightest H 11 
regions el and e2 and the LSB region e3. The blue part of 
the spectrum of region al with hydrogen absorption and 
emission lines is shown in Fig. EJ Table 0] lists the equiva- 
lent widths with their errors of the H7 and H<5 absorption 
lines measured in the wavelength intervals or "windows" 
of Bica & Alloin l|1986f) . The errors include the errors in 
the fitting of line profiles with Gaussians and the noise dis- 
persion in the continuum. A careful placement of the con- 
tinuum level is very important for deriving accurate EWs. 
For this purpose, we choose points in the spectrum free of 
nebular and stellar lines, which were then fitted by cu- 
bic splines. The uncertainties were estimated from several 
different measurements of the equivalent widths of hydro- 
gen absorption lines with independent continuum fittings. 
They are of the same order as the errors in Table 0] 



5.1. Age calibration 

The calibration of the age of stellar populations using 
the equivalent widths of the Ha and H/3 nebular emis- 
sion lines, those of the H7 and H<5 stellar absorption lines 
and the spectral energy distributions is discussed in detail 
in Guseva et al. H2001ll2003all2003b|l . Here we only briefly 
describe these calibrations. 



5.1.1. Balmer emission lines 

The temporal evolution of the Ha and H/? emission line 
equivalent widths depends on the star formation history. 
We consider here the two limiting cases of instantaneous 
burst and continuous star formation models. The equiv- 
alent widths for the instantaneous burst model with a 
heavy element mass fraction Z & /20 are calculated using 
the galactic evolution code PEGASE.2 (Fioc & Rocca- 
Volmerange 1997). The dependence of the Ha emission 
line equivalent width on time is shown in Fig. 6a of Guseva 
et al. l|2003bj) by the thick solid line. These models are ap- 
propriate for regions el and e2 with strong emission lines. 
The equivalent widths of Balmer emission lines in region 
el (EW(Ra) = 998A and EW(H(3) = 166 A) and region 
e2 (EW(Ha) = 872A and EW(R/3) = 134 A) correspond 
to an instantaneous burst age of 4 Myr. 

However, for the LSB regions, models with continu- 
ous star formation are more appropriate. For these mod- 
els we adopt a constant star formation rate (SFR) within 
the time interval from t\ when star formation starts to 
if when it stops. Time is zero now and increases to the 
past. The equivalent widths of hydrogen emission lines 
and SEDs for a set of instantaneous burst models (Fioc & 
Rocca-Volmerange 1997J) are used to calculate the tempo- 
ral evolution of EWs for continuous SF with a constant 
SFR. The temporal dependence of the equivalent widths 
of the Ha emission line is shown in Fig. 6a of Guseva et 
al. I|2003bjl for different U and tf. 

5.1.2. Stellar Balmer absorption lines 

Another way to derive the age of a stellar population is to 
use the relation between the H<5 and H7 absorption line 
equivalent widths and age, derived by Gonzalez Delgado, 
Leitherer & Heckman l)1999|l . Their instantaneous burst 
models predict a steady increase of the equivalent widths 
with age from 1 Myr to 1 Gyr. However, they did not ex- 
tend the calculations for ages ib 1 Gyr when the equivalent 
widths of the absorption lines decrease with age (Bica & 
Alloin 1986). Hence, each value of the hydrogen absorp- 
tion line equivalent width corresponds to two values of 
the age, S3 1 Gyr and <; 1 Gyr. This ambiguity can be 
resolved with the use of other age constraints discussed in 
this paper. 

Furthermore, the models by Gonzalez Delgado et al. 
(1999 ) probably overestimate the equivalent widths of the 
absorption lines at ages ~ 1 Gyr (Guseva et al. [2003b ) . 
Therefore, in the age range from 1 Myr to 16.5 Gyr in- 
stead of the calibration by Gonzalez Delgado et al. {1999 ) 
we use an empirical calibration of the hydrogen absorption 
line equivalent widths versus age by Bica & Alloin (1986). 
This calibration is based on the integrated spectra of 63 
open and globular stellar clusters with known ages, metal- 
licities and reddenings which can be used as templates for 
stellar populations formed in an instantaneous burst. For 
consistency we use the same wavelength intervals or "win- 
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Table 3. Fluxes, equivalent widths of the Ha and H/3 emission lines and the extinction coefficients C(H/3) in the LSB 
regions. 



Telescope 


Region 


Distance a 


Aperture b 


F{Ba) c 


F(H/3) C 








MMT B 


ol 


10.7 


1.5x3.4 


3.53 ±0.14 


1.25 ±0.12 


12.30 ±0.75 


2.72 ±0.24 


0.0 ±0.03 




a2 


29.5 


1.5x13.2 












4m' 


o3 


12.6 


2.0x7.6 


4.88 ±0.39 


1.88 ±0.29 


32.81 ±2.32 


7.72 ±1.01 


0.0 ±0.04 




e3 


18.4 


2.0x2.8 


3.69 ±0.17 


1.19 ±0.14 


109.80 ±1.41 


16.57 ±0.69 


0.04 ±0.06 




o4 


22.6 


2.0x2.8 













Mistance in arcsec from the brightest H II region el. 

b aperture size xxy, where x is the slit width and y the size along the slit in arcsec. 
c in units 10 -16 erg s _1 cm~ 2 . 
d in A. 

e slit orientation with position angle P.A. = 22°. 
f slit orientation with position angle P.A. = 48°. 



Table 4. Equivalent widths of the H7 and US absorption lines in the LSB regions. 



Telescope 


Region 


Distance a 


Aperture b 


EW(HS) C 


EW(H-y) c 


MMT d 


al 


10.7 


1.5x3.4 


-5.53 ±0.22 


-4.92 ±0.23 




a2 


29.5 


1.5x13.2 


-8.89 ±0.53 


-6.07 ±0.41 


4m e 


a3 


12.6 


2.0x7.6 


-8.28 ±0.50 


-5.71 ±0.60 




aA 


22.6 


2.0x2.8 


-9.06 ±0.63 





Mistance in arcsec from the brightest H 11 region el. 

b aperture size xxy, where x is the slit width and y the size along the slit in arcsec. 
c in A. 

d slit orientation with position angle P.A. = 22°. 
e slit orientation with position angle P.A. = 48°. 



dows" for H7 and US flux measurements as Bica & Alloin 
(1986) (A4318-4364A and A4082-4124A, respectively). 

The behaviour of the empirical H<5 absorption line 
equivalent width with the age for an instantaneous burst 
(Bica & Alloin n~98(i|l is shown in Fig. 6b of Guseva et al. 
( 2003b) by the thick solid line. 

The temporal evolution of the H7 and H<5 absorption 
line equivalent widths in the case of continuous SF is 
calculated similarly to that of the Ha and H/3 emission 
line equivalent widths described in the section lB.l.ll More 
specifically, we use the empirical equivalent widths of hy- 
drogen absorption lines (Bica & Alloin 1986.) and SEDs for 
instantaneous bursts (Fioc & Rocca-Volmerange 199TJ) to 
calculate the temporal evolution of EWs in the case of 
continuous SF with constant SFR. The results are shown 
in Fig. 6b of Guseva et al. I|2003b|l for SF with different t\ 
and t[. 



cients obtained from the Ha/H/3 flux ratio in these regions 
are small (Table OJ- However, they are more uncertain as 
compared to the ones in regions el and e2 because of the 
weakness of the emission lines and significant contribution 
of the stellar absorption lines. Ha and H/3 emission lines 
are not detected in the other LSB regions. Therefore, the 
observed SED cannot directly give information on the age, 
but should be used together with the methods discussed 
in Sect . 15 . lTTl and 15". 1 . 21 for simultaneous determination of 
the age and interstellar extinction. 

We used the galactic evolution code PEGASE.2 (Fioc 
& Rocca-Volmerange , 1997) to produce a grid of theoreti- 
cal SEDs for an instantaneous burst of star formation with 
ages ranging between and 10 Gyr, and a heavy element 
mass fraction of Z — Zq/2Q. The SEDs for continuous SF 
in the time interval between t\ ago and tf ago are derived 
by integration of instantaneous burst SEDs. 



5.1.3. Spectral energy distribution 

The shape of the spectrum reflects the properties of the 
stellar population. However, it is also dependent on the 
reddening. A precise determination of the extinction can 
be done only for the two brightest H 11 regions el and 
e2 which possess many strong hydrogen emission lines 
(Table 0. We derived an extinction coefficient C(H/3)=0 
in these regions. In the LSB regions al, a3 and e3, only Ha 
and H/3 emission lines are present. The extinction coeffi- 



5.2. Ages of the stellar populations in the LSB regions 

In this section we derive self-consistently the ages of the 
stellar populations in the LSB regions using: 1) the equiv- 
alent widths of emission lines, 2) the equivalent widths 
of absorption lines, 3) the SEDs and 4) the colours. For 
this we adopt a continuous SF scenario with constant 
or variable SFR. In the latter case we consider a sim- 
plified scenario with two time intervals of SF, which we 
call young and old, with different SFRs. To quantify 
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region aS 

(C(HP)=0.1) 



region aZ 

(C(H0)=O) 




0.0 



4000 5000 6000 

Wavelength (A) 



7000 



4000 5000 6000 7000 

Wavelength (A) 



Fig. 7. Spectra of region a2 on which synthetic SEDs are superposed. Synthetic SEDs are calculated for stellar 
populations forming continuously during one (a) or two (b - g) time intervals. In the case of two intervals the SFR is 
constant within each interval but varies from one interval to another one by a factor b = SFR(young)/SFR(old). The 
spectra in left panel are corrected for interstellar extinction with C(H/3) = 0.1. The extinction coefficient C(H/3) is set 
to in right panel. Time intervals, parameters b and predicted EWs of the Ha emission and US absorption lines are 
indicated (see explanations in Sect. I5.2|i . 



the difference in SFRs we introduce the parameter b = 
SFR(young)/SFR(old). Because of the noisy spectrum of 
region e3 and the significant contamination by nebular 
emission we were not able to measure equivalent widths 
of absorption lines in that region. Therefore, we exclude 
it from the analysis. For the remaining LSB regions, the 
model predictions providing the best fits to the observed 
SEDs are shown in Table For each region we show the 
age range of one or two SF episodes, the parameter b, 
the mass ratio M y /M Q of the young-to-old populations, 
the model equivalent widths of hydrogen lines, the ex- 
tinction coefficients C(H/3) derived from the best match 
between observed and calculated SEDs, and the intrinsic 
and reddened (V — I) colours (see Sect. 15. 2. If) . We use 
the relation E(V - I) = 0.891 C(H/3) (e.g., Aller IT5Mjl 
to correct for reddening. Our young and old stellar popu- 
lations in Table include stars with ages not older than 
250 Myr (log t — 8.4) and not younger than 160 Myr (log 
t = 8.2), respectively. These definitions differ from those 



conventionally used. In fact, the old stellar population in 
our case includes not only several Gyr old stars but also 
intermediate-age stars with age 55 1 Gyr. Negative EWs 
denote absorption lines, positive EWs refer to emission 
lines. Models with highlighted EWs are those in which 
the equivalent widths are not reproduced well despite a 
good fit of the observed SEDs (compare model EWs in 
Table El with observed EWs in Tables Eland EJ). 

We show in Fig. [7| the observed and predicted SEDs 
in the outer LSB region a2 (MMT observations) with no 
emission features. The H7 and US absorption lines in this 
region as well as in the other outer region a4 (4m telescope 
observations) are not contaminated by nebular emission 
from young stellar populations. Hence, the absorption line 
equivalent widths in regions a2 and aA are more accurate 
than in other LSB regions (Table EJ- 

We consider several SF histories and vary extinction 
in the LSB regions to put constraints on the age of their 
stellar populations. Figure[7|shows the spectra of region a2 



Spectroscopic and photometric studies of low-metallicity star-forming dwarf galaxies. III. SBS 1415+437 11 



Table 5. Predicted equivalent widths of emission and absorption lines, extinction coefficients and (V-J) colours from 
models of continuous star formation. 







Age range a 










modeled EW 






modeled (V-I) d 




Region 


young 


old 


b b 


M v /M n 


H5 C 


H 7 C 


H/3 C 


Ra c 




intrinsic 


reddened 






if - 


- ti 


U - 


- U 




















MMT e 


ol 




7.1- 


-7.5 








-6.3 


-5.2 


2.2 


13.8 


0.27 


0.211 


0.452 






7.05 


-7.5 


9.4- 


10 


24 


0.07 


-6.0 


-5.1 


2.5 


13.2 


0.00 


0.456 


0.456 






7.1- 


■7.5 


9.4- 


10 


66 


0.18 


-6.2 


-5.2 


2.2 


12.9 


0.12 


0.335 


0.442 






7.11 


-7.5 


9.4- 


10 


400 


1.00 


-6.2 


-5.2 


1.9 


13.6 


0.22 


0.242 


0.438 




o2 




1.1- 


Q O 

-s.z 








-8.9 


-7.9 


0.6 


3.4 


0.10 


0.350 


0.439 






7.1- 


-8.2 


8.2- 


-9.3 


100 


8.1 


-8.9 


— i .» 


U. ( 


3.8 


0.10 


U.ooo 


U.44 1 






7.1- 


■7.6 


8.9- 


-9.3 


19 


0.5 


-6.8 g 


-5.6 


1.7 


9.7 g 


0.10 


0.342 


0.431 






7.1- 


-8.2 


9.4- 


10 


1000 


19.7 


-8.8 


-7.9 


0.7 


3.9 


0.10 


0.352 


0.441 






7.1- 


-8.2 


8.2- 


-9.3 


7 


0.6 


-9.2 


-7.9 


0.5 


2.7 


0.00 


0.440 


0.440 






7.1- 


-7.6 


8.9- 


-9.3 


9 


0.2 


-7.1 g 


-5.8 


1.4 


7.8 s 


0.00 


0.418 


0.418 






7.1- 


-8.2 


9.4- 


-10 


17 


0.3 


-8.6 


-7.7 


0.6 


3.3 


0.00 


0.441 


0.441 


KPNO 4m' 


a3 




6.7- 


-8.4 








-8.6 


-7.6 


6.6 


36.8 


0.10 


0.355 


0.444 








6.7- 


-8.8 








-9.6 g 


-7.8 


4.3 


22.6 s 


0.00 


0.427 


0.427 






6.7- 


-8.4 


8.4- 


9.3 


7 


1.0 


-8.9 


-7.6 


5.5 


29.3 


0.00 


0.418 


0.418 






6.7- 


-8.4 


9.4- 


10 


20 


0.6 


-8.4 


-7.5 


6.3 


33.0 


0.00 


0.412 


0.412 






6.75 


-8.4 


9.4- 


10 


1000 


32.5 


-8.7 


-7.6 


5.4 


30.5 


0.12 


0.361 


0.468 




aA 




7.1- 


-8.2 








-8.9 


-7.9 


0.6 


3.4 


0.00 


0.350 


0.350 






7.1- 


-8.2 


8.2- 


-9.3 


100 


8.1 


-8.9 


-7.9 


0.7 


3.8 


0.00 


0.358 


0.358 






7.1- 


-7.6 


8.9- 


-9.3 


19 


0.5 


-6.8 s 


-5.6 


1.7 


9.7 s 


0.00 


0.342 


0.342 






7.1- 


-8.2 


9.4- 


-10 


1000 


19.7 


-8.8 


-7.9 


0.7 


3.9 


0.00 


0.352 


0.352 



Muration of SF in log t (t in yr). 
b b = SFR(young)/SFR(old). ' 
C EW in A. 

d intrinsic colour is that for the stellar population in the absence of extinction. (V — J) rc ddcncd = (V — i)i n trinsic + 0.891 (7(14/3). 
°slit orientation with position angle P. A. = 22°. 
f slit orientation with position angle P. A. = 48°. 
g £Ws do not fit observations. 



on which are superposed the synthetic SEDs with different 
SF histories which best fit it. 

We have assumed two extinction coefficients: C(H/3) = 
0.1, derived by Thuan et al. (1999), and corresponding to 
A(V) ~ 0.2, and C{U0) = 0. 

The model SEDs, adopting C(H/3) =0.1 are shown in 
Fig. left panel. The observed properties of region a2 can 
be reproduced by a single young stellar population with 
age between 12 Myr and 160 Myr, if a small amount of 
extinction is assumed (Fig. Up,) . There is no need to invoke 
an older stellar population which, if present, has to be 
much smaller in mass than the young stellar population. 
We note that there is no ongoing SF in this region. The 
most recent SF in region a2 stopped at time 12 Myr ago, 
otherwise its spectrum would have shown a detectable Ha 
emission line (Fig. 6a in Guseva et al. 2003b). 

Next we consider SF scenarios with a varying SFR 
which include older stellar populations with an age of up 
to 2 Gyr (Fig. 03 -Et) and 10 Gyr (Fig. 01). But even in 
these cases, the young population completely dominates 
the light and mass of region a2 (Table 0) . 



An upper limit to the age of region a2 can be obtained 
by assuming no extinction, i.e. C(H/3) = (Fig. right 
panel). Then the observed SED cannot be reproduced by 
a synthetic SED with stars forming between 12 Myr and 
160 Myr ago, as the latter is too blue. However, by vary- 
ing the parameter b, the observed SED can be reproduced 
by SEDs of stellar populations with other SF scenarios. 
In Figs. 0) -01,05 -Ek, we show the best fits labeled by 
the adopted value of b. We exclude the models shown in 
Fig.0: and because they do not fit the observed equiva- 
lent widths of Ha and US. The remaining two models with 
C(H/3) = satisfy all observational constraints. In these 
models, the SFRs for the young population are respec- 
tively 7 times (Fig. Up) and 17 times (Fig.0*) larger than 
SFRs for the old population. These ratios are, however, 
significantly smaller than the corresponding parameters b 
in the case with C(H/3) =0.1. The relative mass fraction 
of the young stellar population is therefore smaller in the 
extinction-free case, being M y /M a — 0.6 (Fig. Up) and 0.3 
(Fig. 00, respectively. Therefore, the presence of a 2 Gyr 
old or even 10 Gyr old population is not excluded in re- 
gion a2, which would however not dominate the optical 
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Fig. 8. (a) V surface brightness distribution along the slit oriented at position angle P. A. = 22°. The origin is set 
at the brightest H II region el. (b) (V — I) colour distribution along the slit with position angle P. A. = 22°. (c) V 
surface brightness distribution along the slit oriented at position angle P. A. = 48°. The origin is set at the brightest 
H II region el. (d) (V — J) colour distribution along the slit with position angle P. A. = 48°. The locations of different 
regions are labeled in (b) and (d). Filled squares show the total modeled colours (stellar plus gaseous emission) for 
region el. Other symbols show the colours predicted for the LSB regions by the models considered (Table EJ). 



emission. If the extinction is non-zero, which is likely the 
case, then there is no need to invoke a significant old pop- 
ulation to explain the observed properties of region a2. It 
is likely that for other regions, al and a3 (Table EJ), the 
extinction is even larger. 

Similar consideration for region a4 shows that the mass 
fraction of the old stellar population is small even in the 
absence of extinction. In any case, the old stellar popu- 
lation, if present does not contribute significantly to the 
luminosity. All the observed properties of this region can 
be reproduced with a young stellar population formed be- 
tween 12 Myr and 160 Myr ago. However, the presence 
of an older stellar population cannot be excluded in that 
case, but one needs to increase the parameter b to match 
the observations (Table [BJ. 

If some extinction is present in region o4, then our age 
upper limits will decrease. There is some hint that the 
extinction may increase with decreasing distance to the 
brightest H n region el. Indeed, by considering models 
which include also a 10 Gyr old population and C(H/3) 
= 0, we find from Table |5] that the closer a region is to 
the brightest part of the galaxy (region al compared to 
region a2, region a3 compared to region o4), the smaller is 



the relative mass M y /M of the young stellar population. 
This is in contrast to the trend found in BCDs where 
the relative mass of the young stellar population decreases 
outwards. Hence, we conclude, that the extinction in SBS 
1415+437 is larger in the inner brighter regions. 



In summary, our spectroscopic analysis of the LSB re- 
gions favors a relatively young luminosity-weighted age 
of the stars populating those regions. A model with a 
constant SFR continuing from 10 Gyr ago until now is 
definitely excluded. An old population > 2 Gyr could be 
present only in models with very specific SF scenarios, 
with a very low SF activity during the first 2-10 Gyr pe- 
riod, a very high star formation rate during the last ~(10- 
200) Myr, and a quiescent period in between. If, however, 
low SF activity has occurred in the period between ^200 
Myr and 2 Gyr ago, then there is no need to invoke a 
stellar population with age Jb 2 Gyr, and all the spectro- 
scopic properties of the LSB regions can be explained by 
only young and intermediate-age stellar populations. 
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Table 6. (V — I) colours of the extended LSB components in some irregular and BCD galaxies. 



No. 


Galaxy 


{V-I) L SB a 


Telescope 13 


12+log(0/H) 


Telescope 


Ref. d 


Ref. c 


1 


NGC 6789 


0.98 


HST 


7.7 or 8.5 


- 


19 


2 


2 


He 2-10 


0.95 


HST 


7.9*" 


4m 


19 


11 


3 


NGC 625 


0.90 


HST 


~8.2* 


1.6m 


19 


16 


4 


Mrk 996 


0.90 


HST 


8.1 


HST 


15 


15 


5 


VII Zw 403 


0.90 


HST 


7.73 


4m 


13 


9 


6 


Mrk 209 


0.88 


2.1m 


7.77 


4m 


19 


8 


7 


UGCA 290 


0.88 


HST 


~7.6 


HST 


19 


1 


8 


NGC 4214 


0.86 


HST 


8.2 


2.1m 


19 


10 


9 


SBS 1437+370 


0.83 


2.1m 


7.93 


4m 


19 


8 


10 


SBS 1319+579 


0.83 


2.1m 


8.11 


4m 


19 


8 


11 


NGC 1705 


0.83 


HST 


8.36 


lm 


19 


14 


12 


NGC 2366 


0.80 


HST 


7.9 


4m 


19 


8 


13 


UGC 4483 


0.70 


HST 


7.58 


4m 


7 


8 


14 


SBS 1249+493 


0.68 


2.1m 


7.72 


4m 


19 


8 


15 


SBS 1331+493 


0.67 


2.1m 


7.77 


4m 


19 


8 


16 


SBS 0940+544 


0.58 


2.1m 


7.50 


Keck 


3 


3 


17 


SBS 1129+576 


0.56 


2.1m 


7.42 


4m 


17 


17 


18 


Tol 65 


0.50 


HST 


7.54 


Keck 


19 


4 


19 


CG 1116+51 


0.50 


2.1m 


7.5 


MMT 


19 


19 


20 


HS 1442+4250 


0.40 


2.1m 


7.63 


4m 


18 


18 


21 


SBS 1415+437 


0.40 


2.1m 


7.61 


4m 


19 


19 


22 


Tol 1214-277 


0.40 


HST 


7.54 


Keck 


19 


4 


23 


Pox 186 


0.40 


HST 


7.73 


MMT 


19 


19 


24 


I Zw 18C 


0.10 


HST 


7.2** 


Keck 


5 


5 


25 


SBS 0335-052 


0.0 


HST 


7.30 


Keck 


12 


6 


26 


I Zw 18 


0.0 


2.1m, HST 


7.2 


Keck 


13 


6 



a colour corrected for reddening in our Galaxy. E(V — i)Gaiaxy are taken from the NED. 

b telescope used to obtain photometry. 

"telescope used to obtain spectroscopic data. 

d reference for the photometry. 

"reference for the chemical abundance. 

*our rough estimate. 

"oxygen abundance of the main component I Zw 18 NW + SE. 

*"oxygen abundance is calculated using the observed data by Kobulnicky et al. (1999) and calibration by Pilyugin (2000). 
References: (1) Crone et al. 2002; (2) Drozdovsky et al. 2001; (3) Guseva et al. 2001; (4) Izotov et al. 2001a; (5) Izotov et al. 
2001b; (6) Izotov et al. 1999; (7) Izotov & Thuan 2002; (8) Izotov & Thuan 1999; (9) Izotov et al. 1997a; (10) Kobulnicky & 
Skillman 1996; (11) Kobulnicky et al. 1999; (12) Papaderos et al. 1998; (13) Papaderos et al. 2002; (14) Storchi-Bergmann et 
al. 1994; (15) Thuan et al. 1996; (16) Vaceli et al. 1997; (17) Guseva et al. 2003a; (18) Guseva et al. 2003b; (19) this paper. 



5.2.1. Age from the colour distributions 

From the HST images, we derived V and / surface bright- 
ness and colour distributions for the regions covered by 
the spectroscopic observations with the MMT (Fig. 
and EJd) and the KPNO 4m telescope (Fig. Et and Hi). 
Note the trend for the (V — I) colour of the LSB com- 
ponent to decrease with increasing distance from region 
el, especially for the KPNO 4m data (Fig. 01). This is 
again suggestive of larger extinction in the brighter re- 
gions. We compare the observed colours with predictions 
from our population synthesis modeling. The results of 
this comparison are shown in Figs. |Sb and IHJi. The pre- 
dicted colours, obtained from convolving the theoretical 
SEDs with the appropriate filter bandpasses, are shown by 
different symbols. The transmission curves for the Johnson 
V and Cousins I bands are taken from Bessell ( 1990). The 
zero points are from Bessell, Castelli & Plez (1998). 



Since the contribution of ionized gas emission to the 
total brightness of region el is significant, the theoretical 
SED for this region has been constructed using a 4 Myr old 
stellar population SED for a heavy element mass fraction 
Z = Zq/20 to which the gaseous continuum SED and 
the observed emission lines were added (see Guseva et al. 
2001). For comparison with the observed colour at P. A. 
= 22°, we reddened the predicted colour of region el by 
an amount corresponding to C(H/3) = 0.11 (Table 3 in 
Thuan et al. 1999), and at P. A. = 48° we adopted C(H/3) 
= (Table 0. The observed colour of region el is very 
blue, (V - I) ~ -0.4 (Fig. and 01), and cannot be 
reproduced by a 4 Myr old stellar population alone ({V 
— I) ~ 0.1). Strong gaseous continuum and line emission 
need to be added (Table ^). 

On the other hand, in all the LSB regions the contri- 
bution of the gaseous emission to the total flux is negligi- 
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Fig. 9. Dependence of the (V — I) colours of the LSB 
component (a) and logarithm of the derived ages (b) on 
the oxygen abundance for the dwarf irregular and BCD 
galaxies from our sample. Dashed lines in (a) show the- 
oretical dependences of the (V-J) colour on the oxygen 
abundance for an instantaneous burst calculated using 
the galactic evolution code PEGASE.2 (Fioc & Rocca- 
Volmerange 1997) in the range of ages (in logarithmic 
scale) log t{ between 7.1 and 10.0 (ti in yr). Solid lines 
show models in which stars are forming continuously at 
a constant SFR between tf =0 and different ti, where log 
ti varies between 7.5 and 10.0. The observed colours are 
corrected only for reddening in our Galaxy. The ages in 
(b) are from Izotov et al. I|2001b}> (I Zw 18), Papaderos 
et al. ltT39~3l) (SBS 0335-05 2), Gus eva et al. lETM!) (SBS 
0940+544), Guseva et al. Il2003bl) (HS 1442+4250) and 
this paper (SBS 1415+437). 

ble. In Figs.lHb and|HJi we show by different symbols the 
colours of the LSB regions with various star formation his- 
tory ( Table EJ|. The total colour (stellar plus gaseous emis- 
sion) of region el is shown by a filled squares. The agree- 
ment between the (V — I) colours and those derived from 
the spectral energy distributions is very good. However, 
some uncertainties in the colours may be introduced by 
the uncertainties in the reddening curves and the ratios of 
total to selective extinction R=A(V)/E(B — V). We use 
R = 3.2 by Aller (1984), based on the reddening curve by 
Seaton (1979). Schlegel, Finkbeiner & Davis (1998) give 
a slightly different R = 3.315. Nevertheless, because the 
extinction in the studied regions is small, the errors intro- 



duced by the use of different reddening curves and R are 
less than 2% in the V band and negligible in the / band. 
The (B — V) colour of ~ 0.2 mag, derived from the SEDs 
of the LSB regions, is also consistent with the observed 
value, derived from the B and V SBPs, which are shown 
in Fig. Hi. 

5.3. Comparison of the LSB component photometric 
properties for galaxies with different metallicities 

Our comprehensive studies of selected galaxies with oxy- 
gen abundances 12 + log(0/H) <^ 7.6 and blue LSB com- 
ponents (SBS 0335-052, Izotov et al. I1997bl Papaderos 
et al. 1X91351 I Zw 18, Izotov et al. 12001b! Papaderos et 
al. TMm SBS 940+5 44, Guseva et al. ETCHl Tol 1 214- 
277, Fricke et al-EDOD Tol 65, P apaderos et al. 113331 SBS 
11 29+576 , Guseva et al. W)o3al HS 1442+4250, Guseva et 
al. I2003bl SBS 1415+437, this paper) have led us to the 
conclusion that these galaxies might be young. This is in 
contrast to the large age estimates for some well-studied 
higher-met allicity irregular and BCD galaxies, such as VII 
Zw 403 and UGCA 290 (Schulte-Ladbeck et al. 1998; 
Crone et al. I2002f) . To investigate this apparent inconsis- 
tency we therefore compare the (V-/) colours of the LSB 
components ((V-I)lsb) of the galaxies from our sample 
with those in galaxies where a large age was derived from 
colour-magnitude diagrams (CMD). 

In Table |HJ we show the oxygen abundances and (V— 
I) colours of the LSB components in several dwarf irreg- 
ular and BCD galaxies. The telescopes used for photo- 
metric and spectroscopic observations are also given in 
the Table. Ground-based photometric data in V and I 
were obtained with the 2.1m KPNO telescope by Y. Izotov 
and R. Green. The HST photometric data were retrieved 
from the archive of the Space Telescope Science Institute 
(STScI) 5 . We included in the sample only galaxies with 
low internal and/or foreground interstellar extinction to 
avoid uncertainties introduced by the correction of the 
(V — I) colour for reddening. Thus, we do not include the 
nearby star-forming galaxy NCG 1569 studied with HST, 
for example. The chemical abundances for the majority 
of galaxies are obtained from spectroscopic observations 
of their H n regions. An exception is I Zw 18C, the faint 
component of the BCD I Zw 18, where no emission lines 
of heavy elements were detected. Therefore, for I Zw 18C 
we adopt the oxygen abundance derived for the bright 
main body of I Zw 18. (V — I) colours of the galaxies 
in Table are corrected for reddening in our Galaxy, with 
E(V — I)Galaxy taken from the NASA/IPAC Extragalactic 
Database (NED). 

The dependence of the (V-I)lsb colours on oxygen 
abundance for 26 selected galaxies is shown in Fig. 
Thin dashed lines represent theoretical dependences for 
an instantaneous burst in the age range from log ti = 7.1 
to log ti = 10.0 (t[ in yr). The models for continuous star 

5 STScI is operated by the AURA, Inc., under NASA con- 
tract NAS5-26555. 
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formation are shown by thick solid lines. These models 
are calculated for a constant SFR which started at time 
t- u with log U between 7.5 and 10.0 and continuing un- 
til now (if = 0). The blueing of the (V — J) colour with 
decreasing oxygen abundance is in agreement with model 
predictions. However, this blueing trend for galaxies with 
12+log(0/H) Ss 7.6 is too steep to be explained only by 
metallicity effects. The blue colours of the LSB compo- 
nents of these low-metallicity BCDs are also not due to 
ionized gas emission, as the latter dominates the outer 
parts of only two BCDs: I Zw 18 (Izotov et al. 2001b; 
Papaderos et al. 2002) and SBS 0335-052 (Izotov et al. 
1997b; Papaderos et al. 1998). In other galaxies, includ- 
ing I Zw 18C, the LSB component emission has mainly a 
stellar origin. The steep trend cannot be explained by red- 
dening effects because the interstellar extinction derived 
from the spectroscopic observations is small for the galax- 
ies shown in Fig.|5Ji. Therefore, it is likely that the blueing 
is mainly due to a change in the age of the stellar popula- 
tions, and that low-metallicity galaxies are younger than 
high-metallicity ones. All objects with 12+log(0/H) Ss 7.6 
and (V-I)lsb Si 0.6 mag are in the range of colours pre- 
dicted for ages Sjl-2 Gyr by instantaneous and continuous 
models of star formation. On the other hand, galaxies with 
(V-I)lsb > 0.7 mag are likely to be older, with ages ~10 
Gyr. Yi (2003) has analysed the uncertainties in the syn- 
thetic integrated (V — I) colours caused by uncertainties in 
the stellar evolutionary models, population synthesis tech- 
niques and stellar atmosphere models. He has shown that, 
despite all the uncertainties, a (V—I) Ss 0.65-0.75 mag can 
be attributed to an intermediate-age stellar population 
with an age not larger than 2 Gyr, in the case of an in- 
stantaneous burst of star formation, in agreement with the 
age derived here from the PEGASE.2 models (dashed lines 
in Fig. EK)- The smaller luminosity- weighted age of galax- 
ies with blue LSB components in FigEK is supported by 
the detailed analysis of the spectroscopic and photometric 
properties of some of these galaxies, as demonstrated in 
this series of papers. 

Izotov & Thuan f!999|) have suggested that the oxy- 
gen abundance may be a good age indicator. That this 
appears to be the case is shown in Fig. Eh> where a clear 
trend of increasing age of the LSB stellar population with 
increasing oxygen abundance is seen for 5 star-forming 
galaxies with available data. The ages for 4 galaxies (I Zw 
18, Izotov et al. 12001 hi SBS 940+54 4, Guseva et al. lSoTffl 
HS 1442+4250, Guseva et al. I2003bl SBS 1415+437, this 
paper) are determined by the four methods described be- 
fore. The age for SBS 0335-052 (Papaderos et al. 11555)) is 
derived using only colours and SEDs. This is because its 
LSB component is embedded into the H II region produced 
by the young central clusters and no absorption lines are 
present in the spectra. 

6. Conclusions 

The results of a detailed photometric and spectroscopic 
study of the metal-deficient blue compact dwarf galaxy 



SBS 1415+437 are presented. Ha images and spectra in 
the optical range have been obtained with the Kitt Peak 
2.1m and 4m telescopes, respectively. A B image has been 
obtained with the 2.2m Calar Alto telescope. These data 
are supplemented by HST/WFPC2 V and I images and 
MMT spectra from Thuan et al. I|1999f) . The main conclu- 
sions of this study can be summarized as follows: 

1. SBS 1415+437 is a nearby (D = 11.4 Mpc) low- 
metallicity cometary BCD with two bright H II regions 
in the SW part of the elongated low-surface-brightness 
(LSB) stellar component. The scale lengths a of the 
LSB component, obtained from surface brightness pro- 
files (SBPs) in V and I are both ~ 0.3 kpc, in excellent 
agreement with those in Thuan et al. I|1999|l . The ob- 
served (V — I) colour of the brightest H II region is 
very blue, ~ -0.4 mag, due to the combined effects 
of a young stellar population and ionized gas emis- 
sion. The colours of the low-surface-brightness com- 
ponent are much redder (~ 0.4-0.5 mag) and roughly 
constant in the outer parts of the galaxy. A deeper 
B band SBP of SBS 1415+437 reveals at intermedi- 
ate and large radii two exponential intensity regimes 
with a scale length of 0.3 and 0.27 kpc, respectively. 
The inner one (4" <, R* <J 13") can be identified with 
the exponential component derived from HST V and 

I data. The outer one, which dominates for R* J> 16", 
is <; 0.3 mag fainter than the inner component. 

2. In the two brightest H n regions we derive oxygen 
abundances of 12 + log(0/H) = 7.61 ± 0.01 and 7.62 
± 0.03 {Z ~ Zq/20). These values agree well with 
previous determinations (Thuan et al. 1999). Other 
heavy element-to-oxygen abundance ratios for these H 

II regions are also in good agreement with mean ratios 
derived from previous studies of BCDs (Thuan et al. 
EEH51 Izotov & Thuan H5J5J) . 

3. The 4 Hc mass fractions Y = 0.246 ± 0.003 and 0.243 
± 0.010, derived for the brightest H n regions of SBS 
1415+437, are in good agreement with previous deter- 
minations for this galaxy and with the primordial 4 He 
mass fraction Y p = (0.244 - 0.245) ± 0.002 (Izotov & 
Thuan IT55B1 Izotov et al. 11559)1 . 

4. We use four different methods and different SF scenar- 
ios to derive the age of the stellar populations in the 
LSB component of the galaxy. For the outer LSB re- 
gions, the equivalent widths of the Hcv and H/3 emission 
lines, the equivalent widths of the H7 and H5 absorp- 
tion lines, the spectral energy distributions and the 
(V — I) colours are reproduced quite well by models in 
which only a young stellar population (t 250 Myr) 
is present. For those regions, only a small extinction 
with C(H/3) in the range 0-0.1 is needed. An older 
stellar population, if present, does not contribute sub- 
stantially to the optical luminosity of those outer re- 
gions. For region al, located closer to the brightest H 11 
region, a larger extinction (C(H/3) = 0.22) is required 
to fit the observational data by the same young stellar 
population. The assumption of C(H/3)=0 in all regions 
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would require an increase of the relative fraction of the 
old stellar population from the outer LSB parts to the 
inner bright H II regions, which would be contrary to 
the trend observed in the majority of BCDs, where the 
relative contribution of old stars increases outwards. 
Assuming no extinction, we find that the upper limit 
to the mass of the old stellar population in SBS 
1415+437, formed between 2.5 Gyr and 10 Gyr, is not 
greater than ~ (1/20 - 1) of the mass of the stellar pop- 
ulation formed during the last ~ 250 Myr. Depending 
on the region considered, this also implies that the star 
formation rate during the most recent star formation 
rate in SBS 1415+437 must be 20 to 1000 times greater 
than the SFR at ages k, 2.5 Gyr. 
5. Wc compare the (V— I) colours of the LSB components 
and the oxygen abundances of a sample of 26 low- 
metallicity dwarf irregular and BCD galaxies, includ- 
ing SBS 1415+437. It is shown that the LSB compo- 
nents are systematically bluer in the lower-metallicity 
galaxies. However, the observed trend is too steep to 
be explained only by metallicity effects. Therefore, 
it is likely that lower-metallicity galaxies have also 
younger populations. In particular, the luminosity- 
weighted ages of galaxies with LSB colours of (V — I) 
Si 0.6 arc probably not greater than 1-2 Gyr. 
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